FRACTAL DIMENSIONS FOR CONTINUOUS TIME RANDOM 

WALK LIMITS 



MARK M. MEERSCHAERT, ERKAN NANE, AND YIMIN XIAO 

Abstract. In a continuous time random walk (CTRW), each random jump follows 
a random waiting time. CTRW scaling limits are time-changed processes that 
model anomalous diffusion. The outer process describes particle jumps, and the 
non-Markovian inner process (or time change) accounts for waiting times between 
jumps. This paper studies fractal properties of the sample functions of a time- 
changed process, and establishes some general results on the Hausdorff and packing 
dimensions of its range and graph. Then those results are applied to CTRW scaling 
limits. 



1. Introduction 

Continuous time random walks have attracted a lot of attention in recent years. 
They provide flexible models for anomalous diffusion phenomena in a wide range of 
scientific areas including physics, finance and hydrology. Consider a random walk 
S(n) = Ji + ■ ■ ■ + J n on M, d , where {J n ,n > 1} model the particle jumps. The con- 
tinuous time random walk (CTRW) imposes a random waiting time between jumps. 
Let T n = W\ + ■ ■ ■ + W n , where {W n , n > 1} are nonnegative random variables. The 
CTRW jumps to location S(n) at time T n . The number of jumps by time t > 
is given by the counting process N t = ma.x{n > : T n < t}, where T = 0. The 
time-changed process S(N t ) represents the location of a random walker at time t > 0. 
A standard assumption in the literature is that {(J n ,Y n ), n > 1} are iid. In recent 
years CTRW with dependent jumps or/and waiting times have also been considered, 
see for example [T21 [251 SI] • 

The scaling limit of a CTRW {S(N t ),t > 0} is a time-changed (or iterated) process 
X = {X(t),t > 0} of the form X(t) = Y(E t ), where the outer process {Y(t),t > 0} 
is the scaling limit of the random walk {S n ,n > 0} and the inner process {E t ,t > 
0} accounts for the random waiting times {W n ,n > 1}. This has been proved by 
Meerschaert and Schefflcr 'I'M, and Becker-Kern, Meerschaert and Scheffler [31 H] 
under the assumption that {( J n , W n ),n > 1} are iid, the jumps { J n , n > 1} belong to 
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Figure 1. Typical sample path of the iterated process X(t) = Y(E t ). 
Here Y(t) is a Brownian motion and E t is the inverse of an 0.8-stable 
subordinator. 

the domain of attraction of an operator stable law and the waiting times {W n , n > 1} 
belong to the strict domain of attraction of a positive stable random variable D with 
index j3 G (0, 1). In this case, the outer process {Y(t),t > 0} is an operator stable 
Levy process with values in M. d and the inner process {(E t ,t > 0} is the inverse of a 
/3-stable subordinator {D(x),x > 0} with D(l) = D. Namely, 

(1.1) E t = inf{x > : D(x) > t}, V t > 0. 

The aforementioned authors further proved that the density function p(t, x) of X(t) = 
Y(E t ) solves fractional partial differential equations; see [TJ [2} |25] and the references 
therein for further information on PDE connections of CTRW limits. 

When the independence assumption on the jumps {J n ,n > 1} is removed, Meer- 
schaert, Nane and Xiao [2B] showed that the outer process Y can be taken as a frac- 
tional Brownian motion, a stable Levy process or a linear fractional stable motion. 
More general inner processes may also be possible if the waiting times are dependent. 

In general, a CTRW limiting process X is non-Markovian, non-Gaussian and sat- 
isfies a form of self-similarity. Figure [1] illustrates a typical trajectory of the time- 
changed process X(t) = Y(E t ), in the case where the outer process Y is a Brownian 
motion. The graph resembles that of a Brownian motion, interrupted by long resting 
periods. This process X is the long-time scaling limit of a CTRW with mean zero, 
finite variance jumps and heavy tailed waiting times in the domain of attraction of a 
/3-stable subordinator, see [29J. 

This paper is concerned with fractal properties of the CTRW limiting process X = 
{X(t),t > 0} defined by X(t) = Y(E t ) for t > 0. In particular we determine the 
Hausdorff and packing dimensions of the range X([0, 1]) = {X(t) : t G [0,1]} and 
the graph GrX([Q, 1]) = {(t,X(t)) : t G [0, 1]}. There has been a large literature on 
sample path and fractal properties of Levy processes [UH H3] , and Gaussian or stable 
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random fields [T§1 HI]. Several methods have been developed for computing the 
Hausdorff dimensions of the range and graph of stochastic processes under "minimal" 
conditions. To give a brief description of the general method, let U = {U(t), t > 0} be 
a stochastic process with values in IR d (for simplicity we assume that the components 
of U(t) are independent). If there exist positive constants C and H G (0, 1) such that 

(1.2) E( sup \U(t + h) -U(t)\ ] < CT H , VTg (0,1), 

V 0<h<T J 

then one can prove 

c\\m„U(\(). m < min ) d. 



1.3) dim H t/([0, 1]) < min <j d, — \ a.s. 



and 

(1.4) dim H Gr?7([0,l]) < minj-^, 1 + (1 - H)dj a.s. 

In the above and sequel, dim H denotes Hausdorff dimension. Moreover, if there exist 
positive constants C and H G (0, 1) such that 

(1.5) &(\U(t) - U(s)\ < \t- s\ H x^j < Cmin{l,x d }, Vs,tG [0,1], 

then equalities hold in both ( II. 3p and (11.41) . The above method can be applied to 
a wide class of stochastic processes, including self-similar processes with stationary 
increments such as stable Levy processes, fractional Brownian motion and the iterated 
Brownian motion ([HIEED])- See [T4] H51 EH] for further information. 

However, the time-changed process X = {Y(E t ),t > 0} considered in this paper 
does not satisfy (11.51) . In fact, if we consider the process X(t) = Y(E t ) in Figure 1, 
where Y is a Brownian motion in W 1 and E t is the inverse of a /3-stable subordinator 
D defined by (II. ip . then the inner process E t remains constant over infinitely many 
intervals, corresponding to the jumps of the stable subordinator D. Hence the graph 
of the process X remains flat over these resting intervals, as evidenced by Figure 1. 
More precisely, it can be proved by using Proposition 2 in Chapter III of Bertoin [S] 
that, for any s < t, F{E S = E t ) > 0. This implies that F{Y(E S ) - Y(E t ) = 0} > 0. 
Hence X does not satisfy (ll.5|) . As we will see from Propositions 13.11 and 13. 6j the 
actual value of dim H GrX([0, 1]) may be strictly smaller than what is suggested by 

d3D. 

Packing dimension was introduced in 1980's by Tricot [12] as a dual concept to 
Hausdorff dimension, and has become a useful tool for analyzing fractal sets and sam- 
ple paths of stochastic processes. It is known that Hausdorff and packing dimensions 
of a set E characterize different geometric aspects of E and many random fractals 
arising in studies of stochastic processes have different Hausdorff and packing dimen- 
sions. See [101 H3] and the references therein for more information and [22l [201 EH] 
for recent development. A fractal set ECtf with the property dim^E = dim p E is 
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usually called a regular fractal. We will see that the range and graph of CTRW limits 
considered in Section 3 are often regular fractals. 

The rest of this paper is organized as follows. In Section 2 we prove under quite 
general conditions that 

(1.6) dim H X([0, 1]) = dim H F([0, 1]) and dim p X([0, 1]) = dim p F([0, 1]), a.s., 

(1.7) dim H GrX([0,l]) = max {l, dim H Z([0, 1])} a.s., 
and 

(1.8) dim p GrX([0, 1]) = max {l, dim p Z([0, 1])} a.s., 

where Z = {Z(x),x > 0} is the R d+1 -valued process defined by Z(x) = (D(x),Y(x)) 
(see (12.121) and (I2.13P below). These results are applied in Section 3 to the scaling 
limits of continuous time random walks. First we consider the uncoupled case, in 
which the iid waiting times {W n ,n > 1} are independent of the iid particle jumps 
{ J n , n > 1}. Then we treat certain coupled examples, where the jump depends on the 
previous waiting time. We also consider triangular array CTRW limits, which lead to 
general inverse subordinators. Finally we examine the case of correlated jumps. In 
all these cases, the outer process Y is either a Levy process or a fractional Brownian 
motion. 

2. General Results 

In this section we prove some general results on the Hausdorff and packing dimen- 
sions of the range and graph of the time-changed process X = {X(t),t > 0} defined 
by X(t) = Y(E t ) for t > 0. We assume that Y = {Y(x), x > 0} is a stochastic process 
with values in R d and E = {E h t > 0} is a process with Eq = and nondecreasing 
continuous sample functions. Both processes Y and E are defined on a probability 
space (Tl, J 7 , P), and they are not necessarily independent. In the following section, 
the results in this section will be applied to CTRW scaling limits. There, E t will be 
taken as the inverse of a strictly increasing subordinator D, defined by (11.11) . For 
a coupled CTRW, where the jump variable depends on the waiting time, the inner 
process E t and the outer process Y(x) in the scaling limit are dependent. 

First we recall briefly the definitions of Hausdorff and packing dimension. More 
detailed information together with their applications to stochastic processes and other 
areas can be found in Falconer [13], Kahane [19], Taylor [ID] and Xiao [13]. For any 
a > 0, the a- dimensional Hausdorff measure of F C M. d is defined by 

(2.1) s a -m{F) = lim inf [ ^(2^)° : F C (J B(x h n), n<e\, 

^ i i=l ' 

where B(x,r) denotes the open ball of radius r centered at x. The sequence of balls 
satisfying the two conditions on the right-hand side of (12.11) is called an e-covering of 

4 



F. It is well-known that s a -m is a metric outer measure and every Borel set in M. d is 
s a -m measurable. The Hausdorff dimension of F is defined by 

dim H F = inf [a > : s a -m(F) = 0} = sup {a > : s a -m(F) = oo}. 

It is easily verified that dim H satisfies the cr-st ability property: For any F n C ~R d , one 
has 

(2.2) dim H ( I J = sup dim H F n . 

V„=i / »>i 

Similarly to (12.11) . the a-dimensional packing measure of F C M. d is defined as 
S a -p(F) = inf|^ S a -P(F n ): Fc|jF n j, 

where s a -P is the set function on subsets of M. d defined by 

s a -P(F) = lim sup< > (2ri) a : B{xi ) r i ) are disjoint, Xj G F, r« < e >. 

The packing dimension of F is defined by dim p F = inf {a > : s a -p(F) = 0}. It 
can be verified that dim p also satisfies the cr-stability property analogous to (12. 2p . 

The packing dimension can also be defined through the upper box-counting dimen- 
sion. For any e > and any bounded set F C R d , let N(F, e) be the smallest number 
of balls of radius e needed to cover F. The upper box-counting dimension of F is 
defined as 

(2.3) dim M F = hm sup . 

£ ^o - log e 



Tricot [12] proved that the packing dimension of F can be obtained from dim M by 
(2.4) dim p F = inf j sup dirn~ M F n : F C (J F n \, 

n n=l ^ 

see also Falconer [T3~| p. 45]. It is well known that for every (bounded) set F C R ', 



(2.5) < dim H F < dim p F < dim M F < d. 

The following theorem determines the Hausdorff and packing dimension of the 
range X([0, 1]) = {X(t) : t G [0, 1]} in terms of the range of Y. 

Theorem 2.1. Let X = {X(t),t > 0} be the iterated process with values in M. d 
defined by X(t) = Y(E t ), where the processes Y and E satisfy the aforementioned 
conditions. If E\ > a.s. and there exist constants c\ and C2 such that for all 
constants < a < oo 

(2.6) dim H F([0, a]) = Ci, dim p F([0, a]) = c 2 a.s., 
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then almost surely 

(2.7) dim H X([0, 1]) = a and dim p X([0, 1]) = c 2 . 

Proof. Since the process t h-> E t is non-decreasing and continuous, the range E([0, 1]) 
is the random interval [0,Ei\. Hence X([0, 1]) = F([0,£i]). 

It follows from the cr-stability of dim H and (12.61) that dim H y([0, oo)) = C\ a.s. Hence 
dim H X([0, 1]) < Ci almost surely. On the other hand, (12. 6p implies 

(2.8) p{dim H F([0,g]) = c 1; Vg G Q + } = 1, 

where Q + denotes the set of positive rational numbers. Since Ei > almost surely, 
we see that there is an event Q' C Q with F(Q') = 1 such that for every w G O' we 
have Ei(u) > and dim H y([0, q],ou) = c\ for all q G Q + . Since for every u £ Q' there 
is a q G Q+ such that < q < Ei(oj), we derive that 

dim H X([0,l],u;) = dim H Y([0, E 1 (co)],u) > dim H F([0, q], u) = c x . 

Combining the upper and lower bounds for dim H X([0, 1]) yields the first equation in 
(12.71) . The proof of the second equation in (12.71) is very similar and is omitted. □ 

Applying Theorem 12. II to the space-time process x i— > (x, Y(x)) with values in IR d+1 , 
one obtains immediately the following corollary. 

Corollary 2.2. Let X = {X(t),t > 0} be the iterated process with values in M. d as 
in Theorem \2.1\ If Ei > a.s. and there exist constants c 3 and c 4 such that for all 
constants < a < oo 

(2.9) dim H GrF([0, a]) = C3 and dim p GrF([0, a]) = C4, a.s., 
then 

(2.10) dim u {(E t ,Y(E t )) : t G [0,1]} = dim H GrF([0, 1]), a.s. 
and 

(2.11) dim p {(E u Y(E t )) : t G [0,1]} = dim p GrF([0, 1]), a.s. 

The random set in the left hand side of (12.101) may be interesting, but it is quite 
different than the graph of X. In order to determine the Hausdorff and packing 
dimension of the graph set of X, we will make use of the M d+1 -valued process Z = 
{Z(x),x > 0} defined on the probability space (Q, J 7 , P) by 

(2.12) Z{x) = (D(x),Y(x)), Vx>0, 
where D = {D(x), x > 0} is defined by 

(2.13) D(x) = inf {t > : E t > x). 

Since t 1— >■ Et is nondecreasing and continuous, it can be verified that the function 
x i—)- D{x) is strictly increasing and right continuous, thus can have at most countably 
many jumps. Moreover, one can verify that D(E t ) > t for all t > and E^^ = x 
for all x > 0. 

6 



Theorem 2.3. Let X = {X(t),t > 0} be the iterated process with values in M. d as 
in Theorem \2.1\ and let Z = {Z(x),x > 0} be the M. d+l -valued process defined by 
\2.12\) and Ii2.13\) . If E\ > a.s. and there exist constants c§ and c$ such that for all 
constants < a < oo 

(2.14) dim H Z([0, a]) = c 5 and dim p Z([0, a]) = c 6 a.s., 
then 

(2.15) dim H GrX([0,l]) = max { 1, dim H Z([0, 1])}, a.s. 
and 

(2.16) dim p GrX([0, 1]) = max {l, dim p Z([0, 1])}, a.s. 

Proof. We only prove (I2.15p . and the proof of (I2.16P is similar. The sample function 
x (—7- D(x) is a.s. strictly increasing and we can write the unit interval [0, 1] in the 
state space of D as 

oo 

(2.17) [0 ) 1] = D([0,£ 1 ))U|J4 

i=i 

where for each % > 1, 1^ is a subintervals on which E t is a constant. Using D we can 
express Ii = [D(xi—),D(xi)), which is the gap corresponding to the jumping site Xi 
of D, except in the case when Xi = E\. In the latter case, Ii = [D(xi— ), 1]. 
Notice that Jj (i > 1) are disjoint intervals and 

Et = E s if and only if s, t G Ii for some % > 1. 

Thus, over each interval Ii, the graph of X is a horizontal line segment. More precisely, 
we can decompose the graph set of X as 

GrX([0,l]) = {(t,Y(E t )):te[0,l]} 

(2 18) °° 

= {(t,Y(E t )) :teD([0, £0)}U \J{(t,Y(E t )) :teh}. 

i=i 

Hence, by the cr-stability of dim H , we have 

(2.19) dim H GrX([0, 1]) = max {l, dim H {(t, Y(E t )) : t e D([0, E 1 ))}}. 

On the other hand, every t £ -D([0, E\)) can be written as t = D(x) for some < x < 
Ei and E t = E D ^ = x, we see that 

(2.20) {(t,Y(Et)) : t e D^Et))} = {(D(x),Y(x)) : x £ [0,^))}, a.s. 
It follows from f l2~T4j) that 

(2.21) p{w:dim H {( J D(x J w),y(z,a;)):xe[0,g])}=C5 J Vg £ Q+} = 1. 
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Combining this with the assumption that E\(u) > almost surely, we can find an 
event Q' 2 ' such that F^fKj') = 1 anc ^ f° r every u G f2 2 ' we derive from (I2.2ip that 

(2.22) dim H { {D{x, co), Y{x, w)):ie[0, £1 (w))} = c B , 

since qi < Ei(u> 2 ) < Q2 f° r some qy, q 2 € Q + , and U C V implies dim H (?7) < dim H (V^). 
Combining (12^U|) and (12T22D yields 

(2.23) dim H {{t,Y{E t )):teD{[0,E 1 ))} = c s , a.s. 

Therefore, (12351) follows from (12TT9D and fl2T23l . □ 

Many self-similar processes Y with stationary increments satisfy (12. 6p and (I2.14p . 
hence Theorems 12.11 and 12.31 have wide applicability. To apply the above theorems 
to CTRW scaling limits, we now take E t to be the inverse of a subordinator D = 
{D(x),x > 0}. We assume that D has no drift, D(0) = and its Laplace transform 
is given by 

E[e- sD(x) ] = e - x ^ D(s \ 

where the Laplace exponent 

POO 

(2.24) *I>d{8)= / (l-e-»)v D (dy). 

Jo 

We also assume that the Levy measure fjj of D satisfies z/d(0, 00) = 00, so that the 
sample function x 1— >• D(x) is a.s. strictly increasing. 

Let Et denote the inverse of D defined by (II. ip . Since the sample function of D is 
strictly increasing, we see that the function t 1— > Et is almost surely continuous and 
nondecreasing. Moreover, P|£?i > 0} = 1. 

Corollary 2.4. Let X = {X(t),t > 0} be the iterated process with values in IR d as 
in Theorem \2.1\ where E t is the inverse (II. ip of a strictly increasing subordinator D 
with -D(O) = 0. Then the conclusions of Theorem \2.1\ Corollary \2.2\ and Theorem 
EJ hold. 

3. Continuous Time Random Walk Limits 

In the following, we compute the Hausdorff and packing dimensions of the range 
and graph of the sample path of scaling limits of continuous time random walks. 

3.1. CTRW with iid jumps: The uncoupled case. Consider a CTRW whose iid 
waiting times {W n , n > 1} belong to the domain of attraction of the positive /3-stable 
random variable D(l), and whose iid jumps { J n , n > 1} belong to the strict domain 
of attraction of the d- dimensional stable random vector Y(l). We assume that {M^n} 
and {J n } are independent; that is, the CTRW is uncoupled. 

It follows from Theorem 4.2 in Meerschaert and Scheffler [29] that the scaling limit 
of this CTRW is a time-changed process X(t) = Y(E t ), where E t is the inverse (11.11) 
of a /3-stable subordinator D. Since D is self-similar with index 1//3, its inverse E is 
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self-similar with index (3. Since Y is independent of E, the CTRW scaling limit X is 
self-similar with index 



Proposition 3.1. Let X = {Y(E t ),t > 0}, where Y = {Y(x) : x > 0} is a stable 
Levy motion of index a 6 (0, 2] with values in M. d and E t is the inverse of a stable 
subordinator of index < (5 < 1, independent ofY. Then 

(3.1) dim H X([0, 1]) = dim p X([0, 1]) = min{d,a}, a.s. 
and 

dim H GrX([0, 1]) = dim p GrX([0, 1]) 

(3.2) J max{l,a} if a < d, 

L + /3(l-I) ifa>d = l, 



a.s. 



Proof. The result ( 13. II) follows from Theorem 12.11 (or Corollary I2.4[) and the known 
results on the Hausdorff and packing dimension of the range of the stable Levy process 
Y. The former is due to Blumenthal and Getoor [6J [T] , and the latter is due to Pruitt 
and Taylor (36]. 

In order to prove ( 13.2[) . we first recall from Pruitt and Taylor [35] their result on 
the Hausdorff dimension of the range of the Levy process Z(x) = (D(x), Y(x)) with 
independent stable components: for any constant a > 0, 

( P ifa</3, 

(3.3) dim H Z([0, a]) = < a if f3 < a < d, a.s. 

( 1 + 0(1-1) if a> d=l, 

Theorem 3.2 in Meerschaert and Xiao [32] (see also Khoshnevisan and Xiao [22] for 
more general results) shows that dim p Z([0, a}) also equals the right hand side of ( 13. 3ft . 
Therefore, (13. 2p follows from the above and Theorem 12.31 □ 

If the CTRW jumps (J n ) have finite second moments, then the limiting process is 
X(t) = B(E t ), where B is a Brownian motion, and Proposition 13.11 with a = 2 gives 

(3.4) dim H GrX([0, 1]) = dim p GrX([0, 1]) = j * + ? )f j > 2 a - s " 

The Hausdorff dimension of the graph of a stable Levy process Y in M. d was deter- 
mined by Blumenthal and Getoor [8j when d — 1 and Y is symmetric, by Jain and 
Pruitt [IT] when Y is transient (i.e. d > a) and by Pruitt and Taylor [35] in general. 
The packing dimension of the graph of Y was determined by Rezakhanlou and Taylor 
I . Combining their results with Corollary 12. 2\ we obtain 

dim H {(£ t ,F(£ t )) : * G [0,1]} = dim p {(£ t , Y(E t )) : t G [0,1]} 
(3-5) J max{l,a} if a < d, 

2 - i if a > d = 1. a ' S ' 

a 
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Clearly, this is different from (13. 2p . Moreover, we notice that the results (13. ip and 
(13. 5 p do not depend on /3, because the set {E t (cj) : t e [0, 1]} IS ct.S. 8u closed interval, 
so the range dimension is the same after the time change. 

3.2. CTRW with iid jumps: The coupled case. In some applications, it is nat- 
ural to consider a coupled CTRW where {(J n ,W n ),n > 1} are iid, but the jump 
J n depends on the preceding waiting time W n . We can also extend the results of 
the last section to certain coupled CTRW limits. In the coupled case, the CTRW 
S(N t ) has scaling limit Y(E t — ) and the so-called oracle CTRW S(N t + 1) has scaling 
limit Y(Et), see [161 HH]- If Y,D are independent, they have a.s. no simultaneous 
jumps, and the two limit processes are the same. The proof of Theorem 12.31 extends 
immediately to the process Y(E t — ) in this case, with the same dimension results. 
This is because the graphs of Z(x) = (D(x),Y(x)) and Z'(x) = (D(x),Y(x—)) have 
the same Hausdorff and packing dimension, as they differ by at most a countable 
number of discrete points. In the following, we discuss examples for Y(E t ), with the 
understanding that the same dimension results hold for Y(E t —). 

The simplest case is W n = J n , in which case X(t) = D(E t ). This process is self- 
similar with index 1, see for example Becker-Kern et al. [3]. It follows from Theorems 
12. 1[ 12.31 and the fact that for any constant a > 0, 

dim H D([0, a]) = dim H {(£>(>), D(x)) : x e [0, a]} = /?, a.s. 

that 

(3.6) dim H X([0, 1]) = dim p X([0, 1]) = 0, a.s. 
and 

(3.7) dim H GrX([0, 1]) = dim p GrX([0, 1]) = 1, a.s. 

We should also mention that by applying the "uniform" Hausdorff and packing di- 
mension results for the /3-stable subordinator D (see Perkins and Taylor [33]), which 
states that almost surely 

dim H D(F) = /3dim H F and dim p D(F) = /3dim p F for all Borel sets FCI, 

we obtain (13. 6p directly by choosing F = [0, Ei). 

Shlesinger et al. [39] consider a CTRW where the waiting times W n > are iid 
with the /3-stable random variable D and E(e _sD ) = e~ sf> and, conditional on W n = t, 
the jump J n is normal with mean zero and variance 2t. Then J n is symmetric stable 
with index a = 2/3. This model was applied to stock market prices by Meerschaert 
and Scalas [2?]. Becker-Kern et al. [3] show that the CTRW limit is X(t) = Y(E t ) 
(t > 0), where Y is a real- valued stable Levy process with index a = 2/3 and E t is the 
inverse of a /3-stable subordinator. Then X(t) is self-similar with index 1/2, the same 
as Brownian motion. However, the Hausdorff dimensions of the range and graph of 
X are completely different than those for Brownian motion. 
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Note that here E t is not independent of Y(t). Theorem 12.11 gives that dim H X([0, 1]) 
= min{ 1,2/3} a.s. To determine the Hausdorff dimension of the graph of X(t), we 
first verify that the Fourier-Laplace transform of (D(1),Y(1)) is 



E ^Y(D-vD(i)^ = E [e" r ' D(1) E(e i « y(1) | J D(l)) 



It follows that the Levy process Z[x) = (D(x),Y(x)) is operator stable [28] with the 
unique exponent 

73- 1 o 

(2f3)-\ 



(3.8) C 



which has eigenvalues /3 _1 and (2/3) ~ 1 . 

By applying Theorem 3.2 from Meerschaert and Xiao [32], we derive that for any 
a > 0, 

(3.9) dim H Z([0, a]) = dim p Z([0, a]) = j ^ f { ^ > J' a.s. 
Consequently, we use Theorem 12.31 to derive 

(3.10) dim H GrX([0, 1]) = dim p GrX([0,l]) =max{l,/3 + i}, a.s., 

which is quite different from the corresponding result (13 .4p in the uncoupled case. 

3.3. CTRW with iid jumps: Triangular array limits. Proposition 13 . 1 1 and ( I3.10p 

rely on the Hausdorff and packing dimension results for sample functions of stable 
or, more generally, operator stable Levy processes. The Hausdorff dimensions and 
potential theoretic properties of general Levy processes have been studied by several 
authors ( [M| |2~T| 123] ) and the packing dimension results have been proved by Khosh- 
nevisan and Xiao [22] and Khoshnevisan, Schilling and Xiao [20]. These results are 
useful for studying fractal properties of the CTRW limits under more general settings 
such as triangular array schemes. 

In particular, for any Levy process Z = {Z(x),x > 0} with values in MP and 
characteristic exponent $ (i.e., E(e*^ ,z ' x ^) = e -21 *^), Corollary 1.8 in [23] shows 
that for any a > 0, 

(3.11) dim H Z([0,a]) = sup \ 7 < p : [ Re ( * ) -ijL < 00 L a.s. 

[ J{teRP: u\\>i} \ L + ll?lr J 

where || • || is the Euclidean norm on MP. On the other hand, Theorem 1.1 in Khosh- 
nevisan and Xiao [22] shows that for any a > 0, 

(3.12) dim p Z([0, a]) = sup S^y > : liminf = J , a.s., 
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where, for all r > 0, W(r) is defined by 
W(r) = 



(More precisely, Theorem 1.1 in [22] is proved for a = 1, but its proof works for 
arbitrary a > 0. Another way to get (I3.12p from Thorem 1.1 in [22] is to use the 
stationarity of increments of Z.) 

By combining (13. lip and (13. 12[) with Theorems 12.11 and 12.31 we extend the results 
in the previous sections to more general time-changed processes. 

Proposition 3.2. Let X = {Y(E t ),t > 0}, where Y = {Y(x) : x > 0} is 
a Levy process with values in M. d and characteristic exponent ip and let E t be the 
inverse of a subordinator D = {D(x),x > 0} with characteristic exponent a. If 
Z = {(D(x), Y(x)), x > 0} is a Levy process in R 1+d and its characteristic exponent 
$ satisfies 
(3.13) 

K~ l Re ( A r^-r ) < Re ( — - ) < KRe 



for all (77, £) G R 1+d with \rj\ + ||£|| large, where K > 1 is a constant. Then almost 
surely, 



< 00 

7 



dim H X([0, 1]) = sup j 7 < : ^ Re « 

and dim H GrX([0, 1]) = max{l,x} almost surely, where 

f , /" / 1 \ d"Qdi 
X = sup < 7 < 1 + a: / Re —r 1 j— < 00 



{ki+ikn>i} vi + ^) + v>(0/(M + ii£ii) 7 

The packing dimensions of X([0, 1]) and GrX([0, 1]) are given as follows, which 
may be different from the Hausdorff dimensions given in Proposition 13.21 

Proposition 3.3. Let X = {Y(E t ),t > 0} be the same as in Proposition Ejj If 
Z = {(D(x), Y(x)), x > 0} is a Levy process in M} +d and its characteristic exponent 
$ satisfies $3.131) . then 

dim p X([0,ll) = sup (7 > : liminf^^ = ol, a.s., 

r->-0+ r 7 J 

where W(r) is defined by 
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W(r) 




1 + a(r)/r) + $i£/r) 



1 



Next we consider a generalized CTRW limit, as in [31], obtained by using a tri- 
angular array scheme. This limit can be applied as a stochastic model for ultraslow 
diffusion (cf. [HI ED- 

At each scale c > we are given iid waiting times (W£) and iid jumps Assume 
the waiting times and jumps form triangular arrays whose row sums converge in 
distribution. More specifically, let S c (n) = J{ + ■ ■ • + J c n and T c (n) =W?+--- + W£, 
we require that S c (cu) =>■ Y{t) and T c (cu) =>■ D(t) as c — )■ oo, where the limits Y(t) 
and D(t) are independent Levy processes. Letting N£ = max{n > : T c {n) < t}, 
the CTRW scaling limit S c (N t c ) => Y(E t ) [311 Theorem 2.1]. 

A power law mixture model for waiting times was proposed in [30J: Take an iid 
sequence of random variables with < Bi < 1 and assume P{W^ C > u\Bi = (3} = 
c" 1 ^ 13 for u > c _1///3 , so that the waiting times are power laws conditional on the 
mixing variables. The waiting time process T c (cu) =>• D(t), which is a subordinator 
with Laplace transform E[e~ sD ^] = e~*^ D ^' ) , where f!2.24p holds. The Levy measure 
of D is given by 



where \x is the distribution of the mixing variable [301 Theorem 3.4 and Remark 5.1]. 



A computation [301 Eq. (3.18)] using - e'^pt-^dt = T(l - /3)s /3 shows that 



Then c 1 A r 4 c E t the inverse subordinator [301 Theorem 3.10]. 

Now we take /i(d/3) = ELi^W 1 " Pk))~%u(W), where < ft <&<■.• < 
/3„ < 1 are constants and 5 a is the unit mass at a. In this case, the subordinator 
is D(t) = Y^k^dkDkit), which is a mixture of independent /3fc-stable subordinators 
Dk{t) (k — 1, • • • , n). See Chechkin et al. [11] for some applications of such CTRW 
and its scaling limit X(t) = Y(E t ). 

In order to apply Propositions 13.21 and 13.31 to establish Hausdorff and packing 
dimension results for the above time-changed process A, we will make use of the 
following technical result. 



(3.14) 




(3.15) 
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Lemma 3.4. Let D(x) = Y^=i dkDk{x), where d k > are constants and D k (x) are 
independent stable subordinators of index fik and < /?i < fa < ■ ■ ■ < fa < 1 . Let 
Y = {Y(x),x > 0} be a strictly stable Levy motion of index a G (0,2] with values in 
IR d . We assume D and Y are independent and let $ be the characteristic exponent 
of the Levy process Z(x) = (D(x),Y(x)). Then for all (77, £) G IR 1+d that satisfies 

H + neii>i 

(3.16, . *.:(, ! ,V K 



where K > 1 is a constant which may depend on n, a, fa, d k for k = 1, 2, • • • , n. 

Proof. For simplicity, we assume that the characteristic exponent of Y is "0(0 = ll£l| a - 
Then the Levy process Z(x) = (D(x),Y(x)) has characteristic exponent 

n 

(3.17) 



k=i 

n 



J2\dkvf k [^s(7T/3 k /2) -ism(nfa/2)\ + 



k=l 

Since fa G (0,1), we have f(r),£) > for all 77, £ G R 1+d . Moreover, < g(rj) < 
Kf(r),£) for some constant K > 0. Hence 

< Re ( ^ — it" ) 



(l + ^ 2 )(l + /(r/,0) " \l + Hv,0J " l + /(»7,0 
From here it is elementary to verify (I3.16p . □ 

By using Lemma T3.4[ Propositions 13.21 and 13.31 we derive the following proposition. 
Since the proof is similar to that of Proposition 4.1 in [32] (see also Proposition 7.7 
in ED), we omit the details. 



Proposition 3.5. Let X = {Y(E t ),t > 0}, where Y = {Y(x) : x > 0} is a strictly 
stable Levy motion of index a G (0,2] with values in M. d and E t is the inverse of a 
subordinator D(t) = ^™ =1 d k D k (t) ; where d k > and Dk(t) are independent stable 
subordinators of index fa an d < fa < fa < ■ ■ ■ fa < 1. Suppose also that E is 
independent ofY. Then 

(3.18) dim H X([0, 1]) = dim p X([0, 1]) = min{d,a}, a.s. 
and 

dim H GrX([0, 1]) = dim p GrX([0, 1]) 

(3.19) J max{l,a} if a < d, 

l + fa(l-l) ifa>d=l, 
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a.s. 



3.4. CTRW with correlated jumps. Now consider an uncoupled CTRW whose 
jumps { J n } form a correlated sequence of random variables, and whose waiting times 
{H^n} are iid and belong to the domain of attraction of a positive /3-stable random 
variable D(l). 

We further assume that { J n } and {W n } are independent. In this case, Meerschaert, 
et al. [26] show that, under certain conditions on the correlation structure of {J n }, 
the CTRW scaling limit is the (if /3)-self-similar process X = {Y(E t ) : t > 0}, where 
Y is a fractional Brownian motion with index H G (0, 1), and E t is the inverse of a 
/9-stable subordinator D which is independent of Y. 

The following proposition determines the Hausdorff and packing dimension of the 
sample path of X. 

Proposition 3.6. Let X = {Y(E t ),t > 0}, where Y is a fractional Brownian mo- 
tion with values in ~R d of index H G (0, 1) and E t is the the inverse of a P- stable 
subordinator D which is independent ofY. Then 

(3.20) dim H X([0, 1]) = dim p X([0, 1]) = min |d, — }, a.s. 
and 

dim H GrX([0, 1]) = dim p GrX([0, 1]) 

(3.21) J i ifl<Hd, 

\ p + (l-Hp)d if I > Hd, a - s - 

Proof. Eq. (I3.20p follows from Theorem 12.11 and the well known result on Hausdorff 
and packing dimension of the range of fractional Brownian motion (see, e.g., Chapter 
18 of [IS])- In order to prove f)3.2ip . by Theorem 12.31 it is sufficient to prove that for 
the R d+1 -valued process Z = {Z(x),x > 0} defined by Z(x) = (D(x),Y(x)),x > 
and for any constant a > 0, we have 

dim H Z([0, a]) = dim p Z([0, a]) 

(3.22) f j. if 1 < Hd, 

\ P + (l- HP)d if 1 > Hd, a ' s ' 

Thanks to ( 12. 5p . we can divide the proof of f)3.22p into proving the upper bound for 
dim p Z([0, a]) and the lower bound for dim H Z([0, a]) separately. These are given as 
Lemmas 13.71 and 13.91 below. □ 

Lemma 3.7. Let the assumptions of Proposition WTb] hold and let a > be a constant. 
Then 

(323) Almp Z( { 0,a])<{l +(1 _ Hg)d f x ll\ a,. 

In order to prove Lemma 13 ,7[ we will make use of the fact that for every e > the 
function Y(x) (0 < x < a) satisfies the uniform Holder condition of order H — e and 
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the following Lemma 13. S\ which is an immediate consequence of Lemma 3.2 in Liu 
and Xiao [21]. It can also be derived from Lemma 6.1 in Pruitt and Taylor [33] . 

Let C7 > be a fixed constant. A collection A (b) of intervals of length b in K is 
called C7-nested if no interval of length b in M can intersect more than C7 intervals 
of A(b). Note that for each integer n > 1, the collection of dyadic intervals I n j = 
[j/2 n , (j + l)/2 n ] is cr-nested with c 7 = 3. 

Lemma 3.8. Let {D(x),x > 0} be a /3 -stable subordinator and let A(b) be a c-j- 
nested family. Denote by M u (b,s) the number of intervals in A(b) which intersect 
D([u, u + s]). Then there exists a positive constant c§ such that for all u > and all 
< b p < s, 

(3.24) e(m u (6,s)) < c 8 sb-P. 
If one takes b = s < 1, then we have 

(3.25) E(M u (a,s)) < c 8 . 
Now we are ready to prove Lemma 13.71 

Proof of Lemma \3. 7\ The proof is based on a moment argument. We divide the 
interval [0, a] into (|_aj + l)2 n dyadic intervals I n j of length 2" n . 

First we construct a covering of the range Z([0, a]) by using balls in IR d+1 of radius 
2~ Hn as follows. Define t n j = j/2 n so that for each I n j = [t n ,j,t n j + 2~ n ], the image 
Y(I n j) is contained in a ball in M. d of radius sup sgJnj ||^(s) — F(t n) j)|| and can be 
covered by at most 

f o at f sn Vsei n . 3 \\Y(s)-Y(t nJ ) 
(3-26) N nJ = eg I 2 _ Hn 



balls of radius 2 Hn . By the self-similarity and stationarity of increments of Y, we 
have 

E(N nJ ) =c 9 2 Hdn E\( sup \\Y(s)-Y(t n>j )\\) d 

(3.27) 



CgE 



( sup \\Y(s 
se[o,i] 



C10 < 00, 



where the last inequality follows from the well known tail probability for the supre- 
mum of Gaussian processes (e.g., Fernique's inequality). 

In order to get a covering for D(I n ^), let r(2~") be the collection of dyadic intervals 
of order n in R + . Let M n j be the number of dyadic intervals in r(2~") which intersect 
D(I n j). Applying ( 13.251) in Lemma [3.81 with b n = s n = 2~ n , we obtain that 

(3.28) E(M nJ ) < c 8 , V 1 < j < ([a\ + 1)2". 

Since 2~ n < 2~ Hn , we see that Z(I n j) = {(D(x),Y(x)) : x 6 I n j} can be covered 
by at most M nJ N nJ balls in R d+1 of radius 2~ Hn . Denote by N(Z([0, a]), 2~ Hn ) the 
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smallest number of balls in R d+1 of radius 2 Hn that cover Z([0, a]), then 

(L«J+1)2" 

N(Z([0,a]),2- Hn ) < Yl MnjNnj. 

3=1 

It follows from (13.271) . (13. 28ft and the independence of Y and D that 
E[iV(Z([0,a]),2- Hre )] < ([a\+l)c 8 c 10 2 n . 

Hence, for any e > 0, 

p{iV(Z([0,a]),2"^ n ) > ([aj + l)c 8 c 10 2"( 1+£ ) } < 2~ ne . 

It follows from the Borel-Cantelli lemma that almost surely 

N(Z([0,a]),2- Hn ) < (|_a| + l)c 8 c 10 2 n ( 1+e > 



K,j = C9 



for all n large enough. This and (I2.3P imply that dim M Z([0, a]) < (l + s)/H a.s. Since 
£ > is arbitrary, we obtain from the above and ( 12. 5 j) that dim p Z([0, a]) < 1/// 
almost surely. 

Next we construct a covering for the range Z([0, a]) by using balls in M. d+1 of radius 
2" n / /3 . Let r(2" n / /3 ) be the collection of intervals in M+ of the form V n k = , 
where fc is an integer. Then the class T(2~ n ^) is 3-nested. Let M' n - be the number 
of intervals in T{2~ n ^) that intersect D(I n J). By Lemma 13.81 with b n = 2~ n ^ and 
s n = 2~ n , we derive E(M' •) < c 8 . Thus D(I n j) can almost surely be covered by ■ 
intervals of length 2~ n ^ from r(2" n / /3 ). 

On the other hand, the image Y(I n j) can be covered by at most 

'sup seIn . \\Y(s)-Y(t n>j )^ d 

2~n/P 

balls of radius 2~ n ^, where t n j = j /2 n , and then similar to ( 13.27ft we derive 
(3.29) HKJ = c 10 2 n{ ?- H)d . 

Denote by N(Z([0, a]), 2~ n /P) the smallest number of balls in lR d+1 of radius 2~ n /' 3 
that cover Z([0, a]), then 

(LoJ+l)2» 

N(Z([0,a}),2-^)< £ KiKr 

By (13.291) and the independence of Y and D we have 

Ek(Z([0,a]),2^)| < (LaJ + l)c 8 c 10 2 n(1+( ?" H)d) . 
Hence, for any e > 0, the Borel-Cantelli Lemma implies that a.s. 

N(Z([0,a]),2- n ^) < (LaJ +l)c 8Cl o2" (1+ ^-^ +e) 
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for all n large enough. This and (12. 3p imply that dim M Z([0, a]) < /3 + (1 — (3H)d + (3e 
almost surely which, in turn, implies dim p Z([0, a]) < j3 + (1 — f3H)d a.s. 
Combining the above we have 



dim p Z([0,a]) < min|— ,/?+(! - /3iJ)d} 



a.s. 



This proves (13T23D . □ 
Lemma 3.9. Under the assumptions of Proposition we have 

(3.30) dim H Z([0.1])>{| +(1 _ WM a., 

Proo/. Since the projection of Z([0, 1]) into R d is F([0, 1]) and dim H F([0, 1]) = ± a.s. 
when 1 < Hd. This implies the first inequality in (I3.30p . 

To prove the inequality in (13.301) for the case 1 > Hd, by Frostman's theorem (cf. 
[T9l p. 133]) along with the inequality 

\\Z(x) - Z(y)\\ > l -[\D{x) - D(y)\ + \\Y(x) - Y(y)\\], 

it is sufficient to prove that for every constant 7 G (0, (3 + (1 — H(3)d), we have 

f 1 f a dxdy 

3.31 E / / ^ < 00. 

Jo Jo [\D(x)-D(y)\ + \\Y(x)-Y(yy" 



Since 1 > Hd, we have /3 + (1 — Hf3)d > d. We only need to verify (13 .3 1[) for every 
7 e (d,p + {l-Hp)d). 

For this purpose, we will make use of the following easily verifiable fact (see, e.g., 
Kahane [T9l p.279]): If S is a standard normal vector in M d , then there is a finite 
constant c\\ > such at for any constants 7 > d and p > 0, 



E 



L(p+I|s| 



< Cn p 



-( 7 -d) 



Fix 1,1/6 [0,a] such that x ^ y. We use Ex to denote the conditional expectation 
;iven the subordinator D, apply the above fact with p = \D(x) — D{y)\ \x — y\~ H and 
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use the self-similarity of D to derive 

E 



1 



[\D(x)-D(y)\ + \\Y(x)-Y(y)\\] 

1 



(3.32) 



\x-y 



< cu\x - y\ 



Ei 



^ 7 E 
1 



| x _ y|-H"(T-rf) 

\D(x) - D(y)\-y- d 



Cl2 



x _ 2/ |Hd+( 7 -d)//3' 



where the last equality follows from the l//3-self-similarity of D and the constant 
C12 = ChE(D(1)~( 7 ~ c( ' ) ) . Recall from Hawkes [T5l Lemma 1] that, as r — > 0+, 



Cl3r M2(l-«) exp 



P(D(1) < r) 
where c 13 = [2tt(1 - £)^/(2(i-fl)] "V^ 



We verify easily Ci 2 < oo. 
It follows from Fubini's theorem and (13.321) that 

dx dy 



E 



(3.33) 



o Jo 



[\Di 



•r 



< C12 



£>(y)l 

dx dy 



\ x _ y \Hd+(y-d)/P 



Y{x) - Y 
< oo, 



the last integral is convergent because Hd + (7 — d)/ (3 < 1. This proves (13. 3 1 f) and 
thus the lemma. □ 

Remark 3.10. Other iterated processes can arise as scaling limits of CTRW with 
dependent and/or heavy-tailed jumps or waiting times. For example, the process Y 
can be taken as a linear fractional stable motion, see [2H]. Our main theorems in 
Section 2 are applicable to these self-similar processes too. However, the problems for 
determining the Hausdorff dimensions of the range and graph sets of the processes Y 
and Z(x) = (D(x),Y(x)) have not been satisfactorily solved, see [381 US] f° r partial 
solutions. 
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